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Two-center two-electron integrals with exponential functions 

Krzysztof Pachuckj*] 

Institute of Theoretical Physics, University of Warsaw, Hoza 69, 00-681 Warsaw, Poland 

Abstract 

We present an efficient approach to evaluate two-center two-electron integrals with exponential 
functions and with an arbitrary polynomial in electron-nucleus and electron-electron distances. We 
show that the master integral with the single negative power of all distances can be obtained from 
the second order differential equation in r, the distance between nuclei. For particular values of 
nonlinear parameters corresponding to the James-Coolidge basis, we find a fully analytic expres- 
sion. For integrals with arbitrary powers of all distances, we construct recursion relations which 
starts from the master integral. The presented approach opens a window for the high precision 
calculations of relativistic effects in diatomic molecules. 
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I. INTRODUCTION 

In order to calculate accurately relativistic effects in atomic or molecular systems, the 
wave function should satisfy cusp conditions. Namely, at the electron-nucleus and electron- 
electron coalescence points the derivative of the wave function is proportional to the wave 
function itself, with coefficients proportional to —Z or 1/2 correspondingly. The widely used 
Gaussian basis set does not satisfy any of these conditions, therefore the numerical conver- 
gence of relativistic effects, is slow or sometimes even does not lead to the right result. It 
has been found recently that inclusion of the single factor e"^''^^ on the top of Gaus- 
sian functions significantly improves convergence of the nonrelativistic energy. Even better 
analytic properties are satisfied by the exponential (Slater) basis set with polynomials of 
electron-nucleus and electron-electron distances. Moreover, the large distance asymptotics 
of exponential functions agrees with that obtained from quasi-classical expansion. However 
the usage of exponential functions in molecular calculations has been limited due to inherent 
difficulties in the accurate and efficient calculations of multi-center integrals. In this work 
we overcome this problem for the simplest case of the two-electron and two-center integral, 
with arbitrary nonlinear parameters and arbitrary powers of electron-nucleus and electron- 
electron distances. We show, that the master integral, with single negative power of all 
distances satisfies a second order differential equation in the nuclear distance r. This equa- 
tion can be solved numerically, or in the case of James-Coolidge basis, is solved analytically 
in terms of Ei- the exponential integral functions. The integrals with higher powers of elec- 
tron distances are obtained by analytic recursion relations which involve the master integral 
and elementary functions. Our approach is based on a set of integration by parts identities, 
and is well established in the area of multi-loop Feynman diagrams . Integration by parts 
identities, similar to those derived here, were recently applied to three-electron one-center 
Hylleraas |3f] and exponentially correlated \^ integrals. 

Let us now define the master two-electron and two-center integral /(r), where r = tab^ 



J An J An Tia ns r2A ru' 

where 1,2 are the positions of the electrons, A,B positions of the nuclei, and the notation 
for nonlinear parameters will be clarified later on. For practical reasons, it is easier at first 
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to consider another integral g defined by 

g{ui) = / / / , (2) 

J An J An J An P23 Psi P12 Pi P2 Ps 

with pi = fi2, p2 = ^2Ai P3 = ^2Bi which is related to / by a Laplace transform, namely 

^/W^= / drf{r)e-'\ (3) 
An r'^ Jo 



2ni J- 



i 00+e 



i 00+e 



The integral (/o = /) with the integer power n > of inter-electronic distance ri2 

d^ri f d^r2 e~"^^^^ g-"2'-i_g g-'!«2f2A g-"'3»'2S ^ 



= / <2-^ (5) 



can be obtained from g by differentiation over Wi, see Eq. fl43|l . A similar situation holds 

for integrals with integer powers i,j, k,l > of all electron-nucleus distances 

/-d^^i . ^3^2 e-"=*''i-^ e-"2^i^ e-'^^^^A g-«'3r-2s ^ 

fn{tj,k,l;r) = / - — / — — ^ — j^. (6) 

J An J An rl/ ry r]/ " 

They can be obtained by further differentiation of g in Eq. ([2]) over nonlinear parameters 

U2,U3,W2,W3. For this we derive recursion relations, which make possible in practice the 

evaluation of integrals with many powers of electron-electron and electron-nucleus distances. 

The function g will be calculated from the pertinent integral in the momentum space, namely 

(7M = G(1, 1,1; 1,1,1), (7) 

where 



2\-m2 



G{mi,m2,m3;m4,m5,mQ) = -^-^ J d^h J d^k2 J d^k-i{kl + u\) ""^{kl + ul 

(kl + uj)-^ {kl2 + wl)-^ {kf, + wl)—^ {kl, + ^3')--(8) 

The topology of these integrals and the notation is presented in Fig. ([T]). 

In Sec. II we will derive a differential equation which is satisfied by the function g. 
This equation can be used to obtain an analytic form of g. However, we find that it is too 
complicated for practical applications. In Sec. Ill we perform an inverse Laplace transform 
to obtain a differential equation for the master integral /(r), which can be conveniently 
solved numerically. In Sec. IV and V, using this differential equation for g, we derive 
recursion relations for evaluation of fn{i,j,k,l;r) in Eq. ([6]). In Sec. VI and VII we work 
out special cases of direct and exchange atomic integrals. In Sec. VIII we consider symmetric 
integrals which are suited for the H2 molecule, and for them we obtain a simple analytic 
form. Finally in Sec. IX we present a short summary. 




FIG. 1: The tetrahedron geometry of the master integral. The double notation is used for vertices, 
namely vertices l(ei) and 0(e2) correspond to the position of the first and the second electron, 
vertices 2 (A) and 3 {B) correspond to the position of nuclei. The nonlinear parameter wi is related 
to ri2 distance and ui = t to r = tab- 

II. DIFFERENTIAL EQUATION 

We use the method of integration by parts identities [2], which is by now standard in 
the analytical calculation of Feynman diagrams. In our case, it amounts to considering the 
following 9 identities in the momentum space representation of the integral G, (i, j = 1, 2, 3) 

d 



= id(i, j) = J d^ki J d^k2 J d^h [k^ {kl + u\ 



{kl + ul)-^ {kl + ul)--^{kl + wlr-^ {kl + wl)--^ {kl, + wl)--^] , (9) 

which are trivially valid, because the integral of the derivative of a function vanishing at 
infinity vanishes. These identities have been first introduced in Ref. 3| for the calculation of 
one-center three-electron Hylleraas integrals. They group naturally into three sets. The first 
set consists of id(l, 1), id(2, 1), and id(3, 1). Other sets are obtained by changing the second 
argument from 1 into 2 or 3. The reduction of the scalar products from the numerator leads 
to identities between functions G of different arguments. Whenever mj = 0, G becomes a 
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known three-body integral, for example 



G(0, 1,1; 1,1,1) ^ 



2wi 



U2 + Us + Wi J \ Wi + W2 + 



2 V ^2 + M3 + 'it'i / 6 



(10) 



where Li is the dilogarithmic function. If we assume all = 1 in Eq. ([8]) and solve 
an arbitrary set of three equations, for example the last one, against three unknowns 
G(l, 1, 2; 1, 1, 1), G(l, 1, 1; 2, 1, 1), ^(1, 1, 1; 1, 2, 1), then the solution for G(l, 1, 1; 2, 1, 1) is 
of the form 

- ^ ^(1, 1, 1; 1, 1, 1) - 2 a ^(1, 1, 1; 2, 1, 1) + P[w^, W2, U2; w^, u^) = 0, (11) 
where a is a polynomial with a tetrahedral symmetry 

222, 222, 222, 222,22/2, 2 2 2 2 2\ 

a = ^2 W3 + U2 % + Ml M3 ^2 + "^1 ""^2 ""^3 + "^1 ""^1 1^1 + — ^2 — M3 — ^2 ~ ^3) 

i22/2|2 2 2 2 2\|22/2|2 2 2 2 2\ /loN 

+U2 W2 {U2 + ^2 - - U3 - ti^i - W3) + U3 ti;3 (M3 + W3 - M2 - «i - w^i - W2), (12) 

and 

P(u;i,ui; W2,U2; ^3,^3) 

= Ml Wi [(Mi + ^2)^ - M3] r(M2 + Wl, Ms, Ml + W2) 
+Ml Wi [{Ui + U'if - wl] T{wi + W73, W2, Ml + M3) 

-[Ui t(7i +u\wl- ul wl + Wl W2 (Mi + M2 ~ ""^3)] r(Mi + M2, Ws, Wi + W2) 

-[uj wf -ulwl+ ul wl + Wl W'i {ul + M3 - wl)] T{ui + M3, W2, Wl + W3) 

+ [m2 {U2 + M^l) (Mi + M3 - wl) - U3 (Mi + M2 - wl)] T{ui + W2, M3, U2 + U^l) 

+ [m3 (^3 + Wl) {u\ + u\- wl) - u\ {u\ + M3 - wl)] T{ui + W^, U2-, M3 + M^l) 

-Wl [W2 {u\ - ul + wl) + W3 {ul +wl- ul)] T{U2 + Us, Wl, W2 + W3) 

-Wl [U2 {ul -wl + ul) + U3 {ul + ul- wl)] T{W2 + W3, Wl, U2 + U3), (13) 



with 



r(„„ „3) = lnK+a.)-ln(a.Y3) . (14) 

(a2 - as) ("2 + as) 



Since G{1, 1, 1; 1, 1, 1) = 5^ and 



G(U.l;2.1.1) = -^^. (15) 
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Eq. (fTT!) takes the form of a differential equation 

0- ^ + ^ ^ 5- + P{Wi, Ui, W2, U2\ W3, M3) = , (16) 

or 

d 

V^T, {\f^g) + P{Wi, Ui] W2, U2; W3, Us) =0. (17) 

OWi 

Analogous differential equations with respect to other parameters can be obtained by ap- 
propriate permutation of arguments, using the tetrahedral symmetry of the function g. The 
solution of this differential equation is presented in the work of Fromm and Hill in in 
the context of the analytic evaluation of the three-electron integral. However, they have 
not used the differential equation approach, but instead performed all integrals directly in 
momentum space. 



III. THE MASTER INTEGRAL 

Let us assume wi = 0. The differential equation in variable t = ui is 

0- + 2 ^ 5' + ^(^' 0; ^3, w^s; ^2, M2) = , (18) 
where a from Eq. (fT2!) is now 

a = [ul-ul + wl- wl) {uj wl - u\ u?^ - {u\ - Ug) {wl - W3) . (19) 
This differential equation takes the form 

{f-j?)g' ^tg = R{t), (20) 



where 



2 ^ {ul -ul^wl- wl) [ul wl - ul wl) 
^ {ul-ul){wl-wl) ' ^ ^ 

P(t,0;M3,W3; ^2,^2) 
{ul - ul) {wl - wl) 

]-\a+ Ri{t) + h+ R2{t) + a- R3{t) + h- Ri{t)\, (22) 



R{t) 



and 



± - ^2 -»3 
a — 2 2 2 2 ' \^'^) 

wi -wi ui- ui 
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Ri{t) 

R2{t) 



U2 



± 



W3 



2 2 2 2 ' 



1 



In 



t + U3 + W3 



t + U-i- W2 V W2 + W3 

m 



t - Ma + W2 



t- U2 + W3 
1 



1 



In 



■"2 + ^3 ^ t + M2 - W3 V W2 + 
t + ^2 + 



t - Ms - W2 Vm2 + M3 + 1^2 + 



t + M3 + U^2 

t-U2- W3 
1 

t + M2 + U^s 



{t + U2 + W2) {t + U3 + W3) {U2 +U3 + W2+ W3) 



(24) 

Inf ^ ' ' "^M , (25) 
Inf ^ ' ' "^M , (26) 

In ^"^ ' '^^ ; "^^^^^ ' -^^^-^Z . -^^ . --.J ^27) 

(t + M3 + U^2) (m2 + tis) iw2 +W3) J ■ 

One performs the inverse Laplace transform of Eq. (!20!) and obtains the differential equation 
for the function / 

r fir) + f{r) - rf{r) + F{r) = , 



{t + U2 + W3){U2 + U3) {W2 + W3) 
t + U3 + W2 



U2+U3+W2 + W3 
{t + U2 + W2) {t + U3 + W3) {U2 +U3 + W2 + W3) 



where 



1 ri cxD+e 

F{r) = —— / dte'''R{t) 

2 IT I J-ioo+e 



a+ Ff (r) + 6+ F2~(r) + a" F^{r) + 6" F4-(r) 

and F~ are the inverse Laplace transform of Ri, 

F^{r) = e"("^-"'2) Ei[-r(u2 + U3)] ± e"^"'^-"^') Ei[-r(w2 + W3)], 
F^{r) = e''("'^-"2)Ei[-r(w2 + w^3)]±e"("^""'^)Ei[-r(M2 + M3)], 

^±(^) = g-r■(^.3+«'2) + W3-U3- W2)iu2 + U3) (^2 + W3) 



(29) 

(30) 
(31) 

(32) 
(33) 

Ei[r(M3 + W2 -M2 - W3)] 

(M2 + U3 + ^i'2 + W3){U2 - U3){W2 - W3) 

+Ei[r{w2 - W3)] + Ei[r{u3 - M2)]} ± e^'("^+'"^) Ei[-r(M2 + % + W2 + ^^3)], (34) 



F±(r) = e-"("2+"'3)|ln 



(U2 + 1^3 - M3 - W^2)(^^2 + U3) {W2 + W3) 



/ , , , w w X Ei[r(M2 + W3-M3-«;2)] 

(U2 + U3 + 1(72 + U^3)(^^2 - U3){W2 " "^3) 

+Ei[r(w3 - W2)] + Ei[r{u2 - M3)]| ± e'^("2+«'3) Ei[-r(M2 + M3 + ^2 + W3)]. (35) 

The solution of the differential equation ( l29l) is 

f{r)=Io{pr) dr'F{r')Ko{pr') + Ko{pr) dr'F{r')Io{pr'), (36) 

where /q and Kq are modified Bessel functions. This is our principal result for the master 
integral /. In this work we do not present any numerical examples for validation of Eq. 
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(136|) . but nevertheless verify, that in the hmit of small r the function /(r) 

/(r) = -rF(0) + 0(r^), 

where 



m 



1 



In 



U2 + U3 



+ 



1 



In 



'W2 + Ws 



(37) 



(38) 



W2 + W:i \U2 + U3 + W2 + W3/ U2+U3 \U2 + U3 + W2 + W3 

coincides with the corresponding helium integral. 

We will show in next sections that all the two-electron two-center integrals can be ex- 
pressed in terms of /, the first derivative /' 



f{r)=p Ii{pr) I dr' F{r')Ko{pr') - Ki{pr) J dr' F{r') Io{pr'] 



(39) 



the exponential integral Ei, and exponential functions. In the derivation of integrals with 
powers of ri2, we will need higher order derivatives Z*-"^ and they can be obtained directly 
from the differential equation fl29l) . 

/(")(r) = (n - 2) ^ /^""'^r) + p^ f''^^\r) - {n - l)^ /^""'^r) - ^ F^''~^\r). (40) 
IV. POWERS OF ri2 



We now pass to the calculation of fnij)) the integral with ^- this we use two 
differential equations with respect to wi and ui 

dg 1 da 



eqi = a —— + - — g + P{ui, wi, U3, W3; W2, M2) = , (41) 
oui 2 oui 

dg 1 da x ^ 

eq2 = ^ ^ ^ o ^ — 9 + P[wi, ui] W2, U2; W3, U3) = 0. (42) 

owi 2 owi 

In the first step the first equation is differentiated n + 1 times with respect to wi and the 
second equation n times, at wi = 0. In the second step, an inverse Laplace transform is 
performed of both equations. According to Eq. ([6]) 



fn{r) = (-1)" 



dwi 



1 



uii=0 2 TT 2 J—ioo+e 



i oo+€ 



dte''' g{t) 



and let us introduce analogous notation 



Un{r) 



Vn{r) 



:-i)" 
:-i)" 
:-ir 



dwi 

Qn 



wi=0 ^.Tfi 



dwi 

Qn 



wi=0 



2ni 



-1 00+e 

i 00+e 



dt e^'' P{t, Wi, U3, W3; W2, U2), 
dt e*'' P(wi, t; W2, U2] W3, %), 



dwi 



01=0 2 7ri J -i 00+e 



/'L 00 + e 
dt e*'' P(W3, U3] W2, U2] Wi,t). 
-i 00+e 



(43) 

(44) 
(45) 
(46) 



In the third step the combination 



d 



— (eqi +req2) +eq2 (47) 

is formed, where all derivatives of /n+i(r) cancel out, and the resulting equation is solved 

against /n+i(r) 

fn+i{r) 



+n {ul - wl) {ul - wl) [r/;_i(r) - (n - 1) /„-i(r)] 
-n {ul + ul + wl + wl) [r /f\(r) + 2 /(^(r)] 
+n[rfi'\{r) + {n + 3)fi%ir)] 
+2 (n - 2) (n - 1) n [r fi%{r) + 2 /fi(r)]|. 



(48) 



This recursion relation allows one to obtain integral with an arbitrary power n of ri2, knowing 
integrals with n — 2, n — 4 and its derivatives with respect to r, for example 
j,^^^^ ^ {ul -ul + wl- wl) {^4 wl - ul wl) f{r) _ a-a+ h' b+ qr f'{r) ^ f"{r) 



p. 



+ [r + 



U2 



+ 



W2 



U2 — ul W2— wl 



-r (U2+W2) 



qr 



+ [r 



U3 



W3 



U2 — ul W2 — wl 



p-r {U3+W3) 



qr 



_ Q-r {U3+W2) _ g-r {U2+W2.) ^+ ^- 

—0 a 



qr 

{U2 - W3) 

b+ 



F+{r) + 



qr 

{W2 + Us) 



2p2 



{W2 - Us) 



F+{r) 



^3^(r) + to)^.(,) 



(49) 



where 



q = {ul-uD {wl-wl). 



(50) 



V. POWERS OF riA,riB,r2A, AND r2B 

Finally we pass to integrals with powers of riA, riB, r2A, and r2B- These are obtained by 
differentiation of /n(r) with respect to corresponding parameters 1^3,1^2,^2, and w^- Let us 
consider differentiation of / with respect to W3. We again use differential equations to derive 
corresponding recursion relations, namely 

d 1 c)(T 

eqi = a- ^-^5- + P{ui, wi, U3, W3; W2, ^2) = , (51) 

oui 2 oui 

dg 1 da . x ^ / x 

^ duT 2 duT ^ ^ ^ ■ (^^^ 



In the first step we differentiate both equations n times with respect to wi and set wi = 0. 
In the second step we perform an inverse Laplace transform. In the third step we form the 
expression 



d / deqi 



+ req3 + eqg 



(53) 



dr V dws 

which cancels out derivatives of /„ with respect to t, and solve the corresponding equation 
against /„(0,0,0, l;r) 

dfn 



/n(0,0,0,l;r) 



1 



2VJr)+rVl(r) + 



{ul -ul + w^- wi) (mI wl - uj wi) 

I / 2 2 4,22,22 r, 2 2\ ri / \ 

-w,{ul~ul) [2/f (r)+r/f (r)] 

9/„_2(r) dfi%{r) 



dws 



+ [n — l)n 



- wl) {ul - w-^) 



dw^ dws 
+r {ul - wl) ws fl_2{r) + 2w3 fi%{r) +rw3 /filr) 



(54) 



In the particular case of n = it takes the form 

Of 



/(0,0,0,l;r) 



dw3 

wi - wi 
1 

2p2 (ty2 _ ^2^ 



a a+wsjul-ul) 

p2 



r fir) 



|a (us - W2) W3 - {u3 + W2) W3 F;^ 



-- [U2 W2 iu2 -Ur^ + W2- W^) - W3 {U2 ^3 - % W^^)] F+ 

+- [u2wl {ul -uj + wl- wj) +W3 {ul w\ - u\w^^\F^\. 
q } 



(55) 



The other single powers of the electron distances can be obtained from the above by ap- 
propriate exchange of U2,U3,W2, and W3. The general recursion can be obtained by further 
differentiation of Eq. fl5^ (after multiplying by the common denominator) with respect to 
W2,U2,W3,U3, or by recursive application of this single differentiation formulae. 

VI. SPECIAL CASE: ATOMIC ORBITALS 



The explicit form of fn{hj, k, I; r) in the general case becomes very lengthy for increasing 
values of i,j, k, and /. Therefore it is worth while to consider special cases which may find 
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practical realization in quantum chemistry codes. When in Eq. ([6]) atomic orbitals are used, 
then two nonlinear parameters in the direct (no exchange) integral are U2 = 102 = 0. We 
thus assume here vanishing of U2 and W2, but allow for an arbitrary polynomial in electron- 
nucleus and electron-electron distances. If we introduce the notation = u, and = w, 



then p = ^/u"^ + w"^, and 



fir) 

and F becomes 
F{r) = 



ri 



^3^2 e""''!^ e-""'^^ 1 1 r 



Att J 4tt riA r2B r2A ru 

poo rr 

loipr) dr' F{r') Koipr') + Koipr) dr' F{r') Io{pr') 
Jr Jn 



(56) 
(57) 



'2m 

„—ur r 



+ 



Ei(— r (u + w)) + Ei(— r u) 
w — u 



2w 



Ei(— r {u + w)) + Ei(— rw) 



2u 



2w 



In 



In 



w + u 
w — u 



w + u 



— Ei(r('u — w)) + Ei(r u) + 2 Ei(— r w) 

— Ei(r(ty — u)) + Ei(r ty) + 2 Ei(— r 



(58) 



All recursion formulae for higher powers of electron distances can be obtained directly from 
the general case considered in the previous section by setting U2 = W2 = 0, and they take 
here a much simpler form. 

VII. SPECIAL CASE: EXCHANGED ATOMIC ORBITALS 



For the exchange integral with atomic orbitals, the relation U'i = W2 = u,U2 = = w 
holds, then p = and using the small x limit of Bessel functions. 



the master integral becomes 
f{r) = 



Io(x) = l+0(x), 
Ko(x) = -(^7 + ln|) +0(x). 

d^ri f d^r2 e""*^''^^^^^^) Q-w(r^B+r2B) 



where 



Atx J Att riAr2A ^15^2^ ri2 

/ dr'Fh') In- = dr' F^-^\r') - , 
Jr r' Jr r' 



F(-^)(r) = / |(e^--e^-)Ei(-2r(M + ^)) 
2 (u^ — w^) I 



(59) 
(60) 

(61) 
(62) 

(63) 
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— e 



+e 



-2rw 



-2ru 



Ei{2r{w - u)) - 2Ei(r {w - u)) + In 
Ei(2 r{u-w))-2 Ei(r {u - w)) + In 



w — u 



w + 


u 


w — 


u 



w + u 



and the superscript ^^^^ denotes not the differentiation, but the integration over r with the 
boundary condition F'^~^\oo) = 0, so dF^~^\r) / dr = F(r). 

The recurrence relations for powers of ri2 is obtained from the differential equation in wi, 
Eq. fj42|) . One divides it by t^, differentiates over wi n-times, performs the inverse Laplace 
transform, and obtains 



fn+iir) 



^ [2n' {u' + w') /„„i(r) - fi'\{r) - (n - 2) (n - if n ^.^ir) 
+Wt'\r)], (64) 



where 



-2ru 



dwi 



wi=0 



(65) 



Wi 



Wi 



-2rw 



r 



Wi — U + W Wi + u + w 



Wi + u — W Wi + u + w 



Wi 



Wi 



2wi 



+ 

r V wi — u + w wi + u — w wi + u + w 
—2 wi Ei(— r {u + w + wi)) . 



(66) 



Similarly differentiation of the master integral with respect to nonlinear parameter, for 
example is obtained from the differential equation in this nonlinear parameter 



/(0,0,0,l;r) 



w 



fir) + 



y2 _ ^2 (y2 



\ (2 Ei(-r iu + w))- e^'"" Ei(-2 riu + w)) 



+e 



-2rw 



Ei(2 r{w-u))-2 Ei(r {w - u)) + In 



w — u 



(67) 



w + u 

The appearance of m — u; in the denominator affects numerical stability of these recursions 
when u ^ w. This problem can be probably solved by employing sufficiently lengthy Taylor 
expansions around u = w, and in this special case / is known analytically, as will be discussed 
in the next section. 
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VIII. SPECIAL CASE: SYMMETRIC 



This is the case when nonhnear parameters are the same for nuclei A and B, namely 
W2 = W3 = w, U2 = U3 = u, and then p = 00. It is the James-Coolidge basis for H2 molecule, 
and was recently used by Sims and Hagstrom for the very accurate calculation of BO 
potential for small nuclear distances. Here we show that all integrals can be expressed in 
terms of Ei and exponential functions. It is convenient in this case to consider a slightly 
different form of the integral, namely 

rfVi f d^r2 e""''!^ g-^'^'is e-'^^^A f,-wr2B ^ 



f{ni,n2,n3,n4,n5;r) 



r2A r2B rl2 



Arc J Ait Tia Tib 
riA - TibT' {r2A - r2BT' {tia + ri^)"^ {r2A + r2B)"^ (68) 



/ for all values of parameters Ui can be obtained from one differential equation in variable 

Ui=t 

0" + 2 5' + ^(^' "3, W3; W2, U2) = (69) 

by the inverse Laplace transform in t and differentiation with respect to Wi, {w2—w^)/2^ {u2 — 
us)/2, {w2 + W3)/2, and {u2+U3)/2, at wi = W2 — W3 = U2—us = 0. This differential equation 
becomes then an algebraic equation, which relates values of / for different arguments, and 
can easily be solved. For example, the master integral is 



/(0,0,0,0,0;r) 



Other examples include 



Ei(-2 r{u + w)) + e-"("+"') (7 + In 



Auw 



2ruw 
u + w 



Ei(-2rw) 



(70) 



/(2,0,0,0,0;r) 



r 



p-r (u+w) ^2 

/(0,0,0,0,0;r) + — 

12uw 

{u + w) e-' -{u-w) e' Ei(-2 r u) 



24 V? w"^ 

-{w-u) e"^ Ei(-2rw;) - (m + w) Ei(-2 r (m + w)) 

+ (« + ^)e-("+-) (7 + ln^'''''^' 
^2 



u + w 



/(0,2,0,0,0;r) = - /(O, 0, 0, 0, 0; r) . 



(71) 
(72) 



Since all other integrals can also be expressed in terms of the Ei and exponential functions, 
matrix elements of the nonrelativistic Hamiltonian can all be obtained analytically. This 
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should allow one to obtain highly accurate nonrelativistic wave functions, and thus precisely 
calculate various relativistic effects to rovibrational energies, shielding and spin-rotational 
constants in the H2 molecule. 



IX. SUMMARY 

We have presented an approach to evaluate two-center two-electron integrals with expo- 
nential functions and with the arbitrary polynomial in electron-nucleus and electron-electron 
distances. All integrals are expressed in terms of the master integral /, the derivative /', 
Ei and exponential functions. The master integral satisfies the second order differential 
equation (!29|) in variable r = tab, and can be accurately solved. This approach certainly 
finds an application in the H2 molecule, for example the present theoretical predictions for 
the dissociation energy 7|, l8| are limited by unknown higher order m corrections and the 
finite nuclear mass effects in the leading relativistic corrections ma^. Both of these cor- 
rections are difficult (if not impossible) to calculate using Gaussian functions. Apart from 
H2, this approach may find applications in arbitrary few electron diatomic molecules. The 
special cases of integrals with direct and exchange atomic functions were considered for this 
purpose. We do not know however, how well relativistic effects can be calculated in this 
aproach. This would require numerical experiments. But the message is that the integrals 
with exponential functions can now be precisely calculated. 
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